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Abstract 

We present a random-matrix analysis of the entangling power of a unitary operator as a function 
of the number of times it is iterated. We consider unitaries belonging to the circular ensembles of 
random matrices (CUE or COE) applied to random (real or complex) non-entangled states. We 
verify numerically that the average entangling power is a monotonic decreasing function of time. 
The same behavior is observed for the "operator entanglement" -an alternative measure of the en- 
tangling strength of a unitary. On the analytical side we calculate the CUE operator entanglement 
and asymptotic values for the entangling power. We also provide a theoretical explanation of the 
time dependence in the CUE cases. 
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I. INTRODUCTION 



In the last years many studies were devoted to the determination of entanglement growth 
laws for bipartite pure states e vol vingfrompro duct states under globally unitary dynamics, 



either with continuous 



naflflflifla 



or discrete time 
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3, 



0, Q|- For not too small systems, and weak couplings between subsystems, the general 
qualitative picture is that of entanglement (subsystem entropy) growing smoothly from zero, 
possibly in a nonmonotonic way, until arriving at an asymptotic regime characterized by 
small oscillations around an equilibrium value. However, when we come to the quantitative 
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level a rich phenomenology is discovered [2i, a, y, i5i, lei, 0, i8i, mi, H2i, iij, uj, Ha m, \m, \m- 

Besides chaos or regularity at the classical level, the choice of parameters like subsystem 
dimensions, coupling strength, initial state, time window, etc., play also important roles in 
determining the law of growth of entropy 191 ]. 

In this paper we concentrate on the regime of very long times, i.e., after the system has 
relaxed to a equilibrium state. More precisely, we are interested in the average value of the 
asymptotic entropy over a suitable distribution of initially nonentangled states. This defines 
the asymptotic entangling power {^(J of the unitary dynamics. 

If the classical dynamics is chaotic in the full phase space, then, according to the Bohigas- 
Giannoni-Schmit conjecture [21I, [3], one should expect that Random Matrix Theory will 
succeed in describing the statistical features of the long time dynamics, in particular, the 
distribution of asymptotic entropies. However, there is a much simpler statistical approach, 
based on the assumption that a typical initial state submitted to a "chaotic" dynamics must 
eventually evolve into a random state, uniformly distributed on the sphere, as far as its 
average properties are concerned. This hypothesis was tested in several finite dimensional 
quantum maps, with a satisfactory quantitative agreement between theory and simulation 
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The purpose of this paper is to compare the predictions of Random Matrix Theory and 
the alluded "Random State Theory" for the average asymptotic entanglement generated by 
a globally unitary map. In Random Matrix Theory the dynamics is explicitly introduced in 
the model: an asymptotic state is generated by the repeated application of a random unitary 
map to a nonentangled initial state [2|. We show that the ensemble of states generated in 
this way does not coincide in general with a uniform distribution on the sphere. 
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Our results are more conveniently stated in a language of operators: the entangling 
power [2^] of U n , where U is a random unitary, decreases (in average) with increasing the 
discrete time n. The statement continues to be true if one substitutes "entangling power" 



by "operator entanglement" 



25 



26 



271 ]. another useful measure of the entangling abilities 
of a unitary (verified numerically, section UTT|) . 

The following section (JTT]) contains the definitions and the exact setting of the prob- 
lem. Sections III 1 1 and lIVI present our numerical and analytical results, respectively. A brief 
discussion of the results is left to section [V] 

II. DEFINITIONS AND SETTING 

We restrict our analysis to the case of bipartite entanglement of pure states in finite 
dimensional Hilbert spaces. As a measure of entanglement, we use the subsystem linear 
entropy. 

Consider a full system divided into two subsystems, A and B. The dimension of the 
full Hilbert space Ti is d = d^dg, with and ds the subsystem dimensions. Let = 
\iPa) <£> \4>b) be a pure separable state of the full system, corresponding to the density 
matrix p = |-0)(-?/>|. In general, after n applications of U, n > 1, the new density matrix 
pi n ) = U n pU r ^ will not correspond to a separable state any more, due to the increasing 
entanglement between the subsystems. This will manifest itself in a growth of the linear 
entropy of the reduced density matrices 

s£\m = l-tr[pW = l-tr[$]\ (1) 

where p^ = tr bP^ and p^ = tr |3- F° r l° n g times, and typical U and \ip), the 
system comes into an equilibrium regime, where the linear entropy shows small fluctuations 
around a stationary average (see, e.g., [12J, 
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23] ), given by 



1 N 

srm)= iim ^E^ n) (i^»- (2) 

n=l 

By doing an additional average on initial product states one arrives at the asymptotic en- 
tangling power of U : 

ep°°(U) = (S?m) m=lM ^ B) ■ (3) 
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It is also useful to consider a time-dependent entangling power, i.e., the initial-state average 
of Eq. ©: 

epM(ZT) = (S^(\m m=lM ^ B) • (4) 

For n = 1 this is just the entangling power of U. 

Concerning the average over product states, we take and \i/jb) to be independent 
random vectors, both of them either real or complex, uniformly drawn from the correspond- 



ing sphere 



29 



30J. In other words, the components of and \4>b) are distributed like the 



columns of a matrix belonging either to the orthogonal group (real case) or unitary group 
(complex case) (Haar measure is assumed in both cases). There are two reasons for these 
choices, (i) They are perhaps the simplest nontrivial cases both from a conceptual point of 
view [29(1 and from the perspective of analytical calculations, (ii) They will allow us to make 
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contact with closely related literature (e.g., Refs. 2I. fl2. 

The problem is how to estimate ep°°(U) for a typical unitary U. For "typical unitary" 
we mean an operator describable (in a statistical sense) by any of the Circular Ensembles of 



Random Matrix Theory (RMT) 3JJ . Accordingly we shall consider that U belongs either to 

the Circular Unitary Ensemble (CUE, unitary group with Haar measure) or to the Circular 

Orthogonal Ensemble (COE), the latter being the appropriate choice for unitaries displaying 

I I 

time reversal symmetry [22J. This leaves us with four cases to analyze: CUE/COE unitaries 
acting on random complex/real states. 

In order to check that our results are not exclusive of the measure chosen for quantifying 
entangling strength [32J, in addition to "entangling power" we also studied the alternative 



measure called operator entanglement 25j] (also known as Schmidt strength 2^]), constructed 
as follows. A bipartite Hilbert space induces a bipartite structure in the space of its linear 
operators, which, equipped with the Hilbert-Schmidt product becomes a bipartite Hilbert 
space itself. Then, operators can be treated as usual vectors, and standard measures for 



entanglement of states can be translated to operators [20| 
linear entropy of the unitary U reads 26] 
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33j. For instance, the 



S L (U) = 1 



1 



E 



ft 1 rf 



(5) 



where the matrix elements of U are related to a product basis, i.e., 



Uk 1 i u k 2 i 2 =a {h\B(h\ U\k 2 )A\i2) 



b ■ 



(6) 
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Of course, by substituting U by U n in the equations above we obtain the operator entan- 
glement as a function of time. 

III. NUMERICAL RESULTS 

We start with a numerical study, emphasizing the most interesting features, but post- 
poning a deeper analysis until section IIVI 

The main ingredients of our simulations are random states (real or complex) and random 



matrices (CUE or COE). They were generated using the same methods as in Ref. 24]. 
Random states are evolved by applying n times a quantum random map. Then, the entropy 
of the final state is calculated and averaged over maps and (if necessary) over states. 

Figure [I] shows the entangling power of a unitary as a function of time, averaged over 
CUE [(a)] or COE [(b), (c)]. The cases (a,b) and (c) correspond, respectively, to complex 
or real initial states. In cases (a) and (c), due to invariance considerations, the average over 
states is redundant and it suffices to consider a single state. For similar reasons, the cases 
n = 1 represent the average linear entropy of standard bipartite pure states, either complex 
[(a)] or real [(b), (c)] (see section [TV]) . 

In the three cases we observe that average entanglement is a decreasing function of the 
number of iterations. This is the opposite to what is observed in weakly coupled maps, i.e., 
entropy increasing from a zero initial value. However, we remark that our purpose is to 
model the equilibrium itself, not the initial phase of relaxation to equilibrium -this would 
require an explicit modeling of the weak coupling, as in Ref. Q • The cases n — 1 and n — > oo 
correspond, respectively, to the predictions of "Random State Theory" and Random Matrix 
Theory. Even though these extreme cases are our main concern, we also analyze the regime 
of intermediate times because it contains valuable information, e.g., about characteristic 
times for the transition between the extremes. 

Evidently the characteristic time for saturation is the Heisenberg time n# = d (in our 
simulations d = 20). In the CUE case the saturation happens abruptly at n = njj. We also 
verified that the operator entanglement behaves in a similar way by plotting the average 
linear entropy of U n clS du function of n (see Fig. [2]). 

Both figures exhibit a very curious characteristic: the asymptotic value for CUE maps 
coincides, within numerical precision, with the n — 1 value for COE. We shall see in the 
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FIG. 1: (a) A fixed initial separable state was chosen arbitrarily and then propagated by U n , with 
U belonging to CUE. Shown is the linear entropy of the evolved states averaged over 10 8 CUE 
matrices, (b) 10 5 COE matrices were used to propagate an ensemble of 10 3 random, complex, 
separable initial states. Shown is the linear entropy averaged over COE matrices and states, (c) 
Idem (a) but for 10 8 COE matrices and one real state. In all cases subsystem dimensions are cIa = 4 
and cLb = 5. Statistical error bars are smaller than symbol diameter. Dotted lines correspond to 
analytical predictions, see section ITVl 

next section that, in the case of the entangling power (Fig. [I]), such a coincidence is indeed 
exact, for all subsystem dimensions cIa and ds- 

IV. ANALYTICAL RESULTS 

The purpose of this section is to explain analytically some of the features present in 
Fig. ([T]) and Fig. ([2]). In some cases we shall be able to understand the global appearance of 
the entangling measures as functions of time, and derive quantitative expressions for some 
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FIG. 2: Shown is the linear entropy of the random operator U n averaged over 10 realizations, 
(a) U belongs to CUE. (b) U belongs to COE. In both cases subsystem dimensions are d A = 4 and 
d B = 5. Statistical error bars are smaller than symbol diameter. Dotted line correspond to the 
analytical prediction, see section IIV1 

limiting values (indicated with dotted lines in Figs. [1] and [2]) . 

The values for the entangling power at n — 1 can already be found in the literature 

ep {1 \U) ia) = d -^ + d ^ + 1 , ( 7 ) 
a + 1 

(1) d 3 -(d A + d B -4)d 2 -[3(d A + d B )-l]d + 2(d A + d B -l) 

6P {U){bUc) = d(d+l)(d + 3) • (8) 



[Subindices (a), (b), (c) refer to each one of the cases depicted in Fig. [TJ] 



above corresponds to the we 



1 known entropy of random complex states 



?he first equality 
1 

301 ] . The second 



result can be found in Ref. {2]. 

Concerning operator entanglement, the case n — 1 was calculated by Zanardi [25| for two 
qudits, i.e., d A = d B - Using techniques to be described below, we obtained the CUE average 
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of Eq. ©: 

d 2 - (d\ + d 2 B ) + 1 
d 2 - 1 



(s L (u)) CVE = , (9) 



thus extending Zanardi's result to arbitrary subsystem dimensions. Inserting dA = 4 and 
dg = 5 in this formula we obtain the value indicated with dotted lines in Fig. [2j 

A. Global features 

All the functions depicted in Figs. [Hand [2] share the property of decreasing in a monotonic 
way and coming to saturation around the Heisenberg time. Or, equivalently, one can say 
that the purity (one minus linear entropy) grows and then saturates. The most surprising 
case is CUE, because of the abrupt saturation at n = n^. This behavior is similar to that 



of the form factor of the Circular Ensembles 22j , defined as 



|trfT| 2 ) = (|t„| 2 ) = J2 (e^-^) , (10) 

a,0=l 

where the average runs over CUE or COE. For CUE the form factor is piecewise linear: 

z 9 \ n if 1 < n < d , 

(\t n \) ={ ~ ~ (11) 

V ' /cuE \d if n>d . 

The explanation for this behavior is as follows. The form factor is a function only of 
the eigenvalues of U n . If n = 1 one has the well known random matrix spectrum which 
shows strong correlations, e.g., level repulsion. For n > 1 the spectrum has been stretched 
and folded n times on the unit circle, and, when n > d the spectrum is almost completely 



uncorrelated [221 ] . Evidently, the same mechanism is responsible for the saturation of the 
entangling measures. 

The similarity between purity and form factor was already noted by Gorin and Seligman 
j^], who considered a continuous time Hamiltonian analogue of the COE case in Fig. [l^c). 
Now we show that such a connection can be established rigorously for CUE maps. Consider 
either the entangling power or the operator entanglement, Eqs. [Hand [5] respectively; insert 
the spectral decomposition for the corresponding unitaries. We recall that eigenvectors and 
eigenvalues are statistically independent in the circular ensembles. In all cases the result 
can be written as follows: 

d 

S(n) = l- C a ^ 7 (e in( ^ + ^-^-^ ) ) . (12) 

0,0,8,1=1 



On the left, S(n) represents any of the average entropies considered. The coefficients C a ^ 
contain the average over eigenvectors and (where applicable) initial states. The time de- 
pendence comes from the average over four eigenphases (p. Due to invariance properties of 
CUE and COE, the averages above do not depend on the particular values of the indices 
a, f3, 5, 7, but only on their being all different, all equal, equal in pairs, etc. Thus, one is left 



with the problem of evaluating a few nontrivial averages 

(exp [in((f) a + - & - 7 )]) , (13) 

(exp [in(2(f>a -<j>s~ 7 )]) , (14) 

(exp [2in(<f> a - &)}) , (15) 

(exp [in(4> a - 4>s)]) • (16) 

For CUE, we can show that all these four averages can be expressed in terms of the basic 
form factors [34( 

(|tn| 2 ) 2 , (U2n| 2 ) , (Un| 2 ) (17) 

[this is immediate for averages ffl5|) and ffl6|) ]. The information we have gathered is enough 
for concluding that in CUE cases one must have 

S(n) = ci + c 2 (|t„| 2 ) + c 3 <J| t 2 n| 2 ) + c 4 (l t n | 2 ) , (18) 

where cj~ are certain time-independent coefficients. 

This result is not unexpected, as the same three basic functions above also appear in the 



CUE average of | t n \ , calculated by Haake some years ago 



221], 



t n | 4 ) = 2 (| t n \ 2 f + (\t 2n \ 2 ) - 2 (| t n \ 2 ) ; (19) 



and (| t n \ A ) is structurally very similar to the entangling measures we are considering: 

d 

(\tX)= E (e m( ^+^-^-^ ) ) . (20) 

a,/3,<5, 7 =l 

Whatever the exact values of Ck in Eq. (T181) . the preceding analysis proves that for CUE 
both entangling power and operator entanglement decay quadratically and then saturate 
abruptly. (Strictly speaking the decay is piecewise quadratic; however this effect is not 



perceptible in our figures, neither is it in a plot of ^| t n \ J versus n 

The possible relationship between the form factor and the entangling measures in the 
COE cases remains a conjecture (Gorin-Seligman's); the required calculations are rather 
more difficult and will not be attempted here. 
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B. Asymptotic values 



As in the preceding subsection, the starting point for the calculations of asymptotic values 
is the general expression f[T2"l) . In the case of the entangling power the coefficients C a ps~/ are 
the result of a double average over eigenvectors |e M ) and initial states [151 ]. 



a 



(e a \i(j)(i/;\es)(ep\tp)(i(}\e 7 ) tr A [tr B (\e a ) (e s \) tr B (|e /3 )(e 7 |)] 



(21) 



The calculation of the asymptotic entangling power requires the time average ([2]) which 
washes out the eigenvalue dependence but enforces the pairing of indices: a = 5 and j3 — 7, 
or a = 7 and (3 = 5. Thus, one arrives at [l5j, |24j 



trA(pl) 2 -E 



{e^)\ 2 [ti A (p a y A )+tr B ( 



a 

PbPb 



where p A and p% stand for the reduced density matrices of the eigenvector \e a ): 



Pa 



tr B \e a )(e a \ , 



p B = tr A |e Q )(e c 



(22) 

(23) 
(24) 



In cases (a) and (c) of Fig. [T]the average over initial states is redundant. It suffices to 
consider just one fixed initial product state. This is due to the invariance of the Haar measure 
with respect to left/right group actions, either for the unitary [(a)] or the orthogonal [(c)] 
group, combined with the fact that local operations do not change the entropy [2J (recall 
that the eigenvectors of CUE and COE are Haar distributed in the unitary and orthogonal 
groups, respectively). So, in cases (a) and (c) we fix the initial state, e.g., = |1)a® |1)b- 
In case (b) we must average a and \iP)b over their respective spheres. In the unitary case 
(a) one has: 



ep°°(I7)(a) = 1- ( E E \Ull,a\ A Ur S ,aUX ta U r , s >, a U; s , >a - 



a r,r f ,s,s f 

E E Wll^Wll^Urs^U^ig^JJr^i.BU*.^ r — 



aj^0 r,r',s,s' 



E E l^ll,o| 2 |^ll,/3| 2 ^rs,o^* s ' iQ f4's',/3^r's,/3 



(25) 



aj^0 r,r',s,s' 

The expression for the orthogonal case (c) is identical to the prece ding one but for U a real 
unitary matrix. The remaining case, (b), will be exhibited in Ref. 341 ] . 



10 



In all cases, the last step is a group average of products of eight matrix elements (not 
always different) belonging to one or two columns, i.e., "one- and two- vector averages of 
monomials of order eight". For such averages we used the powerful diagrammatic method 
devised by Aubert and Lam for the unitary group [35( and adapted by Braun to the orthog- 



onal case [36j. The method is based solely on the unitarity/orthogonality constraint and 
the invariance of the Haar measure under the group actions. It provides explicit expressions 
for some integrals and recurrence relations for others. As the calculations are lengthy but 



otherwise not illuminating, we skip intermediate steps 34| and jump to the final results: 



ep°°(U) {a) = epV{U) mc) , (26) 
d 4 - (d A + d B - 13) d 3 - [12(d A + d B ) - 47] d 2 - 35(d A + d B -l)d 

ep {u)(c) = (dTWTWTWTe) - (27) 

The first line says that the asymptotic average entropy in the unitary ensemble coincides 
with the n — 1 value for COE [see Eq. (jSJ)], for all dimensions d A and d B . This confirms 
the suspicion caused by examining the data in Fig. [TJ However, we have not been able to 
go beyond the mere analytical verification of the conjecture. The deep reasons for such a 
coincidence -if any- remain a mystery. 

The second expression agrees with Gorin-Seligman's calculation j^, who used a different 



method for averaging monomials over the orthogonal group 
For the case (b) we obtained 



e P °°(U\ h) = y , (28) 



with 



X = d 5 + 12d 4 - (4 + 4 - 41) d 3 - (29) 
[ 12(4 + 4) + 2{d A + d B ) - 30 ] d 2 - 
[ 38(4 + d%) + 18 ] d - 16(4 + 4) + 56(cU + d B ) - 40 , 

Y = (d A + l)(d B + l)(d+l)(d + 2)(d + A)(d + Q) . 

One of the advantages of having explicit analytical expressions is that we can now quantify 
the differences between n = 1 and n — > oo. For CUE this represents the difference between 
the predictions of theories based either on random states or on random dynamics. For 
instance, let us consider the scaling with system size, taking for definiteness d A = 2 (which 
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can be thought of as a definition of entangling power according to the multiqubit Meyer- 
Wallach measure {38]). Then, for large ds one has 



ep {1) (U) {&) -ep°°(U) {a) = -L + ... , (30) 

epW(U) {h) -ep°°(U) (h) = 3^ + ... , (31) 

epW([/) (c) -ep 00 ([/) (c) = g^ + ... ■ (32) 
So, the differences are always of second order in the system size. 

V. CONCLUSIONS 

An initially nonentangled state evolving under a globally chaotic dynamics displays 
asymptotically features of randomness. This can be modeled by assuming that the state 
becomes a completely random state, i.e., uniformly distributed on the sphere. Alternatively, 
one can assume that randomness lies in the dynamics, and find out which is the ensemble 
of final states obtained in this way. We showed that both ensembles are different, i.e., the 
dynamics, even if chaotic, does not generate "canonical" random states. When one includes 
in the model the information that states are generated dynamically, the ensemble-average 
entropy decreases due to additional correlations among the state components. This shows 
up as the difference n — 1 vs. n —>■ 00. The effect is relatively small, i.e., second-order in 
system size, but it can be clearly detected in our figures, and might be important for small 
systems. 

A curious byproduct of our studies is the conclusion that the asymptotic entangling 
measures for CUE operators coincide with the respective n = 1 COE cases. Thus, the effect 
of explicitly including the dynamics in the statistical modeling is equivalent to imposing a 
time reversal symmetry. 

Our results contain also a warning against excessively strong interpretations of the 
Bohigas-Giannoni-Schmit conjecture, which associates classical chaos with quantum ran- 
domness. Naively, one may be led to believe that "more chaos always leads to more entan- 
glement". However, if U is classically chaotic, then U n is more chaotic, at least in the sense 
of a higher rate of phase-space mixing. But we have seen here that higher powers of U may 



be less entangling 



39] 
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